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... A recent conjecture regarding the average of the minimum eigenvalue of the reduced density matrix 

00 ' of a random complex state is proved. In fact, the full distribution of the minimum eigenvalue is 

^^ ' derived exactly for both the cases of a random real and a random complex state. Our results are 

^^ , relevant to the entanglement properties of eigenvectors of the orthogonal and unitary ensembles of 

t^ ' random matrix theory and quantum chaotic systems. They also provide a rare exactly solvable case 

(~| ' for the distribution of the minimum of a set of N strongly correlated random variables for all values 

C^ , of N (and not just for large A*'). 
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I. INTRODUCTION 



Entangelement has been studied extensively in the recent past due to its central role in quantum information 

and possible involvement in quantum computation [l|, I4I • It is desirable in many instances to create states of large 

"jrt . entanglement. Measures of entanglement have been studied mostly in the context ofpure bipartite states, where 

the von- Neumann entropy of either subsystem is one of the measures of entanglement [2|. However there exist other 

measures of entanglement as well, e.g. the so called concurrence for two-qubit systems Q. The entanglement of 

C^ random pure quantum states is of interest as they have near maximal entanglement content, especially in the context 

Ch ' of bipartite entanglement [j|. Apart from the issue of bipartite entanglement, statistical properties of such random 

I , states are relevant for quantum chaotic or non-integrable systems. The applicability of random matrix theory and 

'^ ' hence of random states to systems with well-defined chaotic classical limits was pointed out long back @ . They are 

G [ also of relevance to other systems with no apparent classical limit [g, 0, Q • 

In this paper, we focus on a bipartite quantum system. More precisely, we consider a bipartite partition of a 
TV Af -dimensional Hilbert space H^^^^^ as H^^*^) = Ti''^^ ^Hg'K We can assume without loss of generality N < M. 



As an example of such a bipartite system, A may be considered a given subsystem (say a set of spins) and B may 

represent the environment (e.g., a heat bath). Any quantum state \^p) of the composite system can be generally 
f^ ■ written as a linear combination, \ip) — X]i=i Sa=i ^J." \i^) (^ |q;^) where \i^) and \a^) denote two complete basis of 
C ■ ' 7i^ and Ti^ respectively and the coeflficients x^^a's form the entries of a rectangular (iV x M) matrix X . Mutually 

nonexclusive properties of such a state are entanglement, randomness and statisical purity. Such a quantum state lip) 

is: 



,__! ■ • entangled: if not expressible as a direct product of two states belonging to the two subsystems A and B. Only 

[ — I in the special case when the coefficients have the product form, Xi,a = aiba for all i and a, the state \i(j) = \4>^) (E) \<P^) 

^p^ ' can be written as a direct product of two states \4>^) = J2i=i o^iN^) and \4>^) = J2a=i ba\a^) belonging respectively 
J> , to the two subsystems A and B. In this case, the composite state \ip) is fully unentangled. But otherwise, it is 
r ""j ■ generically entangled. 

^ '■ • random: if the coefficients Xi^a are random variables drawn from an underlyting probability distribution. The 

5t 1 simplest and the most common random state corresponds to choosing Xi^a's as independent and identically distributed 
Gaussian variables, real or complex. 

• pure: if the density matrix of the composite system is simply given by, p = |V^)(V'| with the constraint Tr[p] = 1, 
or equivalently (■0|'0) = 1- 

Given a random, pure and generically entangled composite state, important informations on the results of the 
measurement of any observable on the subsystem A can be derived from the reduced density matrix pA = Tr^ [p] , 
obtained upon tracing out the environmental degrees of freedom (i.e., those of subsystem B). It is easy to show 
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(see section II for details) that for a random pure state, pA = XX'^ is an N x N square matrix where X is the 
N xM rectangular coefficient matrix. The A^ unordered eigenvalues Ai, A2, • . . , A^r of pA carry important informations 
regarding the degree of the entanglement in the subsystem A. Given that the entries Xi^a of the coefficient matrix X 
are independent Gaussian variables (real or complex), the eigenvalues A^'s of the matrix pA = XX^ are also random 
variables and their joint probability density function (jpdf ) is known [9|, llO| 



^ ^, 



\j=l / j=l j<k 

Here /3 = 1,2 corresponding to the real and complex entries of A and Bm,n is the normalization constant that is 
known explicitly [lO|. Several spectral properties associated with the jpdf in Eq. ([T]), in particular for the complex 
(3 = 2 case, have been studied extensively in the literature, for instance see the book [ll| and references therein. 

In principle, all informations about the spectral properties of the subsystem A^ including its degree of entanglement, 
are encoded in the jpdf ([T]). For example, one useful measure of entanglement is the von Neumann entropy S = 
— ^j^-^ Ai In(Ai) which is a random variable. The average entropy {S) (where the average is performed with the 
measure in Eq. ([1])) was computed for /? = 2 by Page [l2| and was found to be {S) ~ ln(A'^) — -^ for large 
1 << N < M . Noting that In(A^) is the maximal possible value of entropy of the subsystem A, it follows that in 
the limit when M » N, the average entropy, and hence the average entanglement, of a random pure state is near 
maximal. Later, the same result was shown to hold for the /3 = 1 case [13| . 

While the average entropy is a useful measure of entanglement, it is not the unique one. In fact, important 
informations regarding the nature of entanglement of a random pure state can also be obtained (see Section II for a 
detailed discussion) by studying the probability distributions of the extreme eigenvalues Amax = max(Ai, A2, . . . , Xn) 
and Amin = min(Ai, A2, . . . , Aat). In particular, the probability distribution of the minimum eigenvalue Amin provides, 
in addition to the nature of the entanglement, an important information about the degree to which the effective 
dimension of the Hilbert space of the subsystem A can be reduced. 

In fact, the avera ge v alue (Amin) (with respect to the measure in Eq. ((T|)) of the minimum eigenvalue was studied 
recently by Znidaric [Ij] for the case N = M and based on the exact (Amin) for small values of N, Znidaric conjectured 
that (Amin) — l/N^ for all N for the complex case (/? = 2). The purpose of this paper is to provide exact results for 
the full probability distribution of Amin for all N (for the case when N = M), both for the complex {(3 = 2) and the 
real (/3 = 1) cases. A byproduct of our general results is the proof of Znidaric's conjecture for /? = 2. Our results are 
summarized as follows. Let PN{x)dx denote the probability that x < Amin < a; + dx, i.e., Pn{x) is the probability 
density function (pdf) of Amin- We show that 

• Complex case (/? = 2): 

Pn{x) = TV (TV^ _ 1) (1 _ Nx)'^'-^ 9(1 - Nx) (2) 

where Q{x) is the standard Heaviside function, Q{x) = 1 for x > and 0(a;) = for a; < 0. The k-th moment 
HkiN) = (A^i„) is given by 

_ r(fc + i)r(iv^) 

In particular, for fc = 1, we get pi{N) = 1/N'^ thus proving the recent conjecture in [14| . 

• Real case (/3 = 1): the result for the real case turns out to be a bit more complicated. For the pdf of Amin we 
get 

P.ix) = A. x--l^ (1 - ^.)(-^--)/^ 2F, r ^, ^, ^-1±1-^^ -\^1±\ ,0<x<l/N (4) 

and Pn{x) = for x > 1/N . The constant An is given by 

jvr(iv)r(ivV2) 

N 2^-ir(iV/2)r((7V2 + iV-2)/2)' ^' 

and 2i^i(a, /?, 7, z) is the standard Hypergeomctric function defined as [15| 

a/3 , a(a + l)/3(/? + 1) z2 a(a + l)(a + 2)/3(/3 + l)(/3 + 2) z^ 



The moments fik{N) = (A^j^) are also computed exactly and are given in Eq. (ISS)) . In particular, the average value 
(fe = 1) decays for large N as 



A^i(^)-]^ 



(7) 



where the prefactor c has a nontrivial value 



c = 2 1 



= 0.688641 • ■ 



(8) 



The paper is organized as follows. In Section II, we provide a general introduction to the random pure states of a 
bipartite system and recapitulate some general facts leading to the jpdf H]). Section II and III provide the detailed 
calculations of the distribution of the minimum eigenvalue for the complex and the real cases respectively. Finally we 
conclude in Section IV with a summary and open questions. Some details of the calculations are presented in the two 
appendices. 

II. A RANDOM PURE STATE OF A BIPARTITE SYSTEM 

In this section we recall some general facts about a random pure (RP) state of a bipartite system, its entanglement 
properties and the associated random matrix ensemble. As mentioned in the introduction, let us consider a composite 
bipartite system A i^ B composed of two smaller subsystems A and B, whose respective Hilbert spaces 7i^ and 

Tig have dimensions N and M. The Hilbert space of the composite system H^^*^-' — Ti)^ H^ is thus NM- 
dimensional. Without loss of generality we will assume that N < M. Let {|«^)} and {[a^)} represent two complete 
basis states for A and B respectively. Then, any arbitrary state \tp) of the composite system can be most generally 
written as a linear combination 



N M 
i—l Q — 1 



a 



(9) 



where the coefficients Xi^a's form the entries of a rectangular {N x M) matrix X = [xi^a]- 

Now, the state \tp} is a statistically pure state of the composite system if the density matrix of the composite system 
is given by 



p = m{^\- 



(10) 



Note that had the composite system been in a statistically mixed state, its density matrix would have been of the 
form 



k 



(11) 



where \ipk)^s are the pure states of the composite system and < p^ < 1 denotes the probability that the composite 
system is in the fc-th pure state, with ^j, p^ = 1. In this paper, we will restrict ourselves to the case when the composite 
system is in a pure state denoted by lip). Then its density matrix in Eq. (jlOp . upon using the decomposition in Eq. 
([9]), can be expressed as 



i,a j,l3 



(12) 



where the Roman indices i and j run from 1 to A^ and the Greek indices a and (3 run from 1 to M . We also assume 
that the pure state jV') is normalized to unity so that Tr[p] — 1. Hence the coefficients XiaS must be such that 
Tr[p] = 1. 

Given the density matrix of the pure composite state in Eq. (|12p . one can then compute the reduced density matrix 
of, say, the subsystem A by tracing over the states of the subsystem B 



M 



PA 



TrB[p]-5](«^|p|«^ 



(13) 



a=l 



Using the expression in Eq. ((T2| one gets 

N M N 

PA^T. E^^"4o N^)0-'l = E W^,\^^){A (14) 

i~j — l a—1 ^,J — 1 

where Wij's are the entries of the N x N square matrix W = XX^ . In a similar way, one can express the reduced 
density matrix pg = Tr^[p] of the subsystem B in terms of the square M x M dimensional matrix W — X'^X. 

Let Ai, A2, . . . , Xn denote the N eigenvalues oi W = XX"^ . Note that these eigenvalues are nonnegative, Ai > for 
alH = 1, 2, . . . , iV. Now the matrix W' — X^X has M > N eigenvalues. It is easy to prove that M — TV of them are 
identically and N nonzero eigenvalues of W^ are the same as those of W. Thus, in this diagonal representation, one 
can express pA as 

N 

p^ = ^A.|Af)(Af| (15) 

where |A^)'s are the eigenvectors oi W — XX^ . A similar representation holds for pB- It then follows that one can 
represent the original composite state lip) in this diagonal representation as 

N 

l^>=EVA;|Af)®|Af) (16) 

where |A^) and |A^) represent the normalized eigenvectors (corresponding to nonzero eigenvalues) of M^ = XX'^ and 
W' = X^'X respectively. This spectral decomposition in Eq. (fT6|) is known as the Schimdt decomposition. The 
normalization condition (iplip) — 1, or equivalently Tr[p] — 1, imposes a constraint on the eigenvalues, X]i=i ^i ~ ^■ 

Note that while each individual state \Xf) (X) \Xf) in the Schimdt decomposition in Eq. (I16|) is unentangled, their 
linear combination \ip), in general, is entangled. This simply means that the composite state \ip) can not, in general, 
be written as a direct product \ijj) = \(j)^) (g) |(/>^) of two states of the respective subsystems. The spectral properties 
of the matrix W, i.e., the knowledge of the eigenvalues Ai, A2, . . . , Aat, in association with the Schimdt decomposition 
in Eq. (|16p . provide useful information about how entangled a pure state is. For example, as mentioned in the 
introduction, one useful measure of the entanglement is the von Neumann entropy, S — — X]i=i ^i lii(-^i)- 

In addition, the two extreme eigevalues, the largest Amax — max(Ai, A2, . . . , Aat) and the smallest Amin = 
min(Ai, A2, . . . , Ajv) also provide useful information about the entanglement. Note that due to the constraint 
^1=1 ^i = 1 ^ii'i the fact that all eigenvalues are nonnegative, it follows that 1/N < Amax < 1 and < 
Amin < 1/A^. Consider, for instance, the following limiting situations. Suppose that the largest eigenvalue 
Amax = max(Ai, A2, . . . , Ajv) takes its maximum allowed value 1. Then due to the constraint X]z=i -^i — 1 ^^"^ 
the fact that Ai > for all i, it follows that all the rest {N — 1) eigenvalues must be identically 0. In that case, it 
follows from Eq. ([16]) that jV') is fully unentangled. On the other hand, if Amax = 1/^ (i-e., it takes its lowest allowed 
value), it follows that all the eigenvalues must have the same value, Ai = 1/A^ for all i, again due to the constraint 
Si=i Xi = 1. In this case, one can show that the pure state {tp) is maximally entangled, as this state maximizes the 
von Neumann entropy S — In(A^). 

In this paper, we will focus on the smallest eigenvalue < Amin < 1/A^. As in the case of the largest eigenvalue 
above, let us consider the two limiting situations. When Amin takes its maximal allowed value Amin = 1/-^, it follows 
again from the constraint J2i=i ^i ~ ^ that all the eigenvalues must have the same value Ai — 1/N. This will thus 
make the state \ip) maximally entangled. In the opposite case, when Amin = takes its smallest allowed value, while it 
does not provide any information on the entanglement of the state IV'), one sees from the Schmidt decomposition that 
the dimension of the effective Hilbert space of the subsystem A gets reduced from A^ to iV — 1. Indeed, if Amin is very 
close to zero, one can effectively ignore the term containing Amin in Eq. I|16p and thus achieve a reduced Hilbert space, 
a process called 'dimensional reduction' that is often used in the compression of large data structures in computer 
vision [1^, [13, [i3l • Thus the knowledge of Amin and in particular its proximity to its upper and lower limits provide 
informations on both the entanglement phenomenon as well as on the efficiency of the dimensional reduction process. 

So far, our discussion is valid for an arbitrary pure state in Eq. ^ with any fixed coefficient matrix X = [xi^a]- 
Now, such a pure state will be called a random pure state if the coefRcients a^i^a's are random variables, drawn 
from an underlying probability distribution. In particular, we will consider the case when the elements of X are 
independent and identically distributed random variables, real or complex, drawn from a Gaussian distribution: 



Prob[X] ex exp -|Tr(A:tX) 



, where the Dyson index j3 — \,2 corresponds respectively to the real and complex 



X matrices. The product W = XX^ is called the random Wishart matrix [T^. The joint distribution of the A^ 
nonnegative eigenvalues of W is known [20| 

•8=1 j<k 

Note however, that in case of a random pure state |-0) in Eq. ©J the eigenvalues of the matrix W — XX"^ are not 
quite the same as that of the Wishart matrix, due to the additional constraint that Tr[/9] — Tr[VF] = 1. Thus, the 
eigenvalues of W that appear in the Schimdt decomposition in Eq. P^ . arc distributed according to the Wishart law 
in Eq. (|17p . but in addition have to satisfy the constraint X]j=i ^i = 1- This constraint can be explicitly incorporated 
by multiplying a delta function (5(X]i=i ^i ~ 1) to the Wishart measure in Eq. pT|) . With this additional delta function 
multiplying the Wishart measure, the exponential term in Eq. (|17p just becomes a constant and can be absorbed 
into the overall normalization constant and one arrives at the jpdf of the eigenvalues of W mentioned in Eq ([1]) in the 
introduction. 

Given the jpdf ([T]), we are interested here in the distribution of the minimum eigenvalue Amin- Let Qn,m{x) — 
Prob[A,„in > x] be the cumulative distribution of Amin- The pdf of Amin is simply obtained by taking the derivative, 
Pn,m{x) = —dQ]S!,M{x)/dx. Since the event Amin > x necessarily implies that all the eigenvalues \i > x (for all 
i = 1,2,..., N), it follows, upon using the explicit jpdf H]), that Qn,m{x) is precisely given by the multiple integral 
(with N < M) 

poo poo f N \ ^ li. 

The real technical challenge is to evaluate this multiple integral. In the next two sections, we show how to compute 
this integral exactly respectively for /? = 2 and (3—1, for all M — N, i.e., when the Hilbert spaces of the two 
subsystems have equal dimensions. In this case, i.e., when M = N, we will denote, for simplicity of notations, 
Qn,n{x) = Qn{x) for the cumulative distribution of the minimum eigenvalue and the corresponding density by 

Pn.n{x) = Pn{x) ^ -dQN[x)/dx. 

III. A COMPLEX RANDOM VECTOR 

This section is devoted to finding exactly the distribution of the minimum eigenvalue Amin or the minimum Schmidt 
coefficient for random complex states. Let 

Qn{x) = Prob [Amin >x\= Prob [Ai > 2:, A2 > 2:, . . . , Aat > a;] . (19) 

Therefore 



Qn{x)^Bn,n •••/ 5[Y^K~l\\[{\j-\kfJ[d\, 
"'^ "'^ \i=i J j<k 1=1 



(20) 



3< 

An evaluation of this multiple integral proceeds by introducing an auxiliary one defined by 



poo poo / N \ N 

/ ■■•/ 5 EA.-t ^(^^-^^•)'^'^^- 



nx,t)= .../ 5\i^K-t\[[ix,-Xkr[[d\, (21) 

"'^ "'^ \i=i / j<k i=\ 

so that Qn{x) = Bn,n I{x, 1). Consider the following Laplace transform of I{x, t): 

/•OO /"OO ^ 

I{x,t)e-'*dt^ I •••/ e-'^"-^^^l[{X,^\k)^l[dX,. (22) 

A linear shift and scaling Zi — s{Xi — x) results in 



" •'^ •'■^^ j<fe i=i 



S JQ JQ 



I{x,t)e-'*dt^—j^ 1^ ■I e-^^-^'^l[{z,-Zkyl[dz,. (23) 

j<k i=l 



Thus the dependence on s and x just factors out of the integral. The integral happens to be one of the Selberg 
integrals which can be evaluated explicitly [21[ and this gives 






(24) 



An inverse Laplace transform yields 



nf-"o' r(7 + 2)r(j + 1) ^2 , 

I{x,t)^ -^^■'-° ':^^' {t-Nxf -^ eit-Nx). (25) 



Using the known normalization constant [lO| 

r(7v2) 



we finally arrive at 



Bn.n = — Kf — i (26) 

nf=o r(iv-j)r(iv-j + i) 



Qn{x) = Prob [Amin > a;] = Bn,n I{x, 1) =- (1 - Nx)^ "^9(1- Nx) . (27) 



Subsequently, the pdf is given by 



dQN{x) _ ^ji^j2 i\/i AT \JV^-2, 



Pn{x) = ^^f^-^ = N{N' - 1)(1 - NxY^ ^' 9(1 - Nx). (28) 

ax 

A plot of this pdf can be found in Fig. 1 for A^ = 4. Thus Pn{x) in a; G [0, 1/A^] has the limiting behavior 

Pn{x) -* N{N^ -I) as a; ^ 

= N{N^ -l){l~-Nx)^^-^ as x^l/N (29) 

Note that in the regime where x << 1/N, the pdf in Eq. (P5|) becomes exponential, Pn{x) « N{N'^ — 1) exp[— A(A^ — 
l)a:]. Let us also note that the distribution of the smallest eigenvalue in Eq. (|?7|) is identical to that of the smallest 
intensity component of a complex random state derived recently |22l |. provided one replaces A^^ (in the exponent in 
Eq. (EZD) by A^. 



Moments of Amin: From the explicit expression of the pdf in Eq. ([28)1 one can easily compute all the moments of 
Amin- For the fc-th moment we get 

M.(A^) = (AL„> = l^ x^Pn{x) dx = ^^\^^^l^fj^ - (30) 



In particular, for fe = 1, we obtain for all A^ 



^ll{N) = (A„,in) = ^, (31) 



thus proving the recent conjecture by Znidaric |14l | based on evaluations for small A^. Putting fc = 2 in Eq. (|30|) . we 
get the second moment ^2 = ivrnv^+TT" Thus the variance is given by 

^'-^^2{N)-[,,{N)Y = j^,[j^- (32) 

IV. A REAL RANDOM VECTOR 

While complex random vectors are "generic", real vectors are important as well. For instance in the case when 
the system has a time-reversal symmetry or any anti-unitary symmetry the eigenfunctions can be in general chosen 
to be real and the relevant ensembles are the "orthogonal" ones (such as the Gaussian orthogonal ensemble and 
the circular orthogonal ensemble), wherein general orthogonal transformations leave the ensemble invariant [2ll |23||. 
The entanglement properties of real and complex random states may, in general, differ. For instance for so called 



"single-particle" states or one-magnon states, real states have lower entanglement measured in terms of two-spin 
entanglement content than the case of the complex states [2J]. In general, much less is known for random real states 
than the complex ones, although for instance several many-body Hamiltonians (say of spins) have natural time-reversal 
symmetry. In this section the distribution of the minimum eigenvalue of the real case is calculated exactly. 
The jpdf of the eigenvalues Xi in this case (we again restrict ourselves to the case M — N) is 



/ AT \ N 

Piv(Ai, • • • , \n) = Cn,nS ^ a. - 1 m |Aj - Afcl n -7^, 



(33) 



where Cn n is the normalization constant and is known to be llC 



^N,N 



,Ar N-l 



Cn]n - ^^/2r(jv2/2) n r [^-^) r [^-^j ■ (34) 

The cumulative distribution of the smallest eigenvalue, Qn{x) = Prob[Ainin > x], is given by 



poo noo f N \ AT 



QNix)^CN,N I •••/ 5| >^A,-1| II |A,-Afc|||^=dA,;. (35) 

To evaluate this multiple integral, we proceed, as in the previous section, by defining an auxiliary integral J{x, t) 
as 

Ji^^t)^ I ■I S{yX,-t]Y\\^,-Xk\T]^dK, (36) 



I.OO poo / N \ ^1 



so that Qn{x) — Cn.n J{x, 1). 

Taking the Laplace transform of Eq. (|36p leads to 



-1 poo /"OO ^ -| 

J{a:,t)e-^'dt=-—^ ■■■ e'^^^-^y^ Y[\y, - y,\l[-=dy„ (37) 



y^S) J2sx J2sx .^^, j^;^ Vy- 

where the scaled variable yi — 2s\i. We next use a result due to Edelman \2^ for the Wishart orthogonal ensemble 
whose jpdf is given by 



where the normalization constant a^ ^ is 



i=l V ^* j<fe 



«^.^ = 2^V2r(7V2/2)- (39) 

For such an ensemble Edelman [2^ showed that the distribution of the smallest eigenvalue Q^{z) — Prob [j/min > z] 
is given explicitly by 

JV 



Nr{N) /■~e-^''/2 /iV-l 1 y' 

Myi ■ ■ ■ uyN = -^ 

i=l 

where U{a, b; z) is the confluent hypergeometric function [26| of the second kind that satisfies the differential equation 



«»,.,^f.^^fp"<..^...,n*-».^,.^?§^r9^K^.4^f)* <«) 



with the boundary conditions 






Uia,b,0)^ J^}] ^\ C/(a,5,z^oo)-0. (42) 

i (1 -|- a — o) 



Working back we therefore obtain 



J{x,t)e-''^dt 



(2s)^V2 



Nr{N) 



2N-i/2Y{N/2)aN,N 



-'e-^y/\ fN-l 1 



2sx 



Vv 



y 



(43) 



To make further progress, it turns out to be easier to work with the probabihty density function rather than the 
cumulative distribution Qn{x), 



dQ^ix) dJ{x,l) 

Fn[X) = ; = -Cn,N- 



dx 
Taking the derivative of Eq. (|43)) with respect to x leads to 



dx 



dJ{x,t) t 



dx 



e-^'dt = b 



1 e 



-Ns 



N- 



where 



bN 



/^g{N^-l)/2 



NT{N) 



U 



N -1 1 



2 ' 2'""^ 



r(iV/2)aAr,Ar2(^'+2^-2)/2 



The task then is to find the Laplace inverse: 

_ dJ{x,t) _ 6jv 1 
dx \fx ^ 



-Nsx 



JN^-l)/2 



u 



N ~1 1 



-, sx 



First, an application of the convolution theorem leads to 



C-' 



-Ns 



JN^-l)/2 



N'-l 



{t-Nx)^'^'~^^^^e{t-Nx). 



(44) 



(45) 



(46) 



(47) 



(48) 



Second, using an integral representation of the hypergeometric function U(a, b, z) [26J | namely 

1 



U{a, b, z) 



r(a) Jo 



e-''H''-'{l + t) 



b-a-l 



dt 



one obtains the following inverse: 



C7^ U 



N -1 1 



, , sx 

2 ' 2' 



„3/2 



r(V 



— i(^-3)/2(^^^)-(JV+2)/2 



(49) 



(50) 



Using the two inverses in Eqs. (l48l) and ((50| and the convolution theorem, we get upon simplifying 



C: 



-Nsx 



,(7V2-l)/2 



U 



N -1 1 



, sx 



r3/2 



t-Ni 



j.-{N-l)/2 



p |^ Af^+Af-2 ^ 



, iV^ + JV-4 

(t-Nx) 5 2^1 



t'^ {x + 1') ^ {t-t' - Nx) ^ dt' (51) 

A^ + 2 A^-l N'^ + N-2 t~Nx" 



Here 2Fi{a,b,c,z) is the standard hypergeometric function [2g|, and the integral can be found in [T^. Using this 
along with Eqs. (|44l47p and substituting t = 1, we finally get the p.d.f. of the minimum eigenvalue Amin as 



Pj,{x) = Aj, x-^/2 (1 - iVa;)(^'+^-4)/2 ,F, f 



N + 2 iV-1 N^ + N -2 1-Nx 



2 ' 2 ' 
and Pn{x) = for x > 1/N . The constant An is given by 

NV{N)T{N^I2) 



A 



N 



■)N-l 



T{N/2)T{{N^ + N -2)12) 



,Q<x<l/N (52) 



(53) 



This solves exactly for the distribution of the minimum eigenvalue of the reduced density matrices of bipartite 
random real states when the dimensions of the subspaces are equal. In the simplest possible case of real states of two 
qubits, N — 2, the distribution is simply 



P2ix) 



\-2x 



V^ 



< a; < 1/2; ^2(2;) = 0, x > 1/2. 



(54) 



This follows from Eq. ([5^ as 2^1(2, 1/2, 2, x) — l/-y/l — x. Alternatively it almost immediately follows from the jpdf 
in Eq. (I33p as there are only two eigenvalues that sum to unity in this case, and the distribution of the one which is 
less than one-half is precisely P^ix)- In Fig. 1, we plot the pdf Pn{x) of A,nin for A'' = 4, both for the complex case 
given in Eq. ([28| and the real case given in Eq. ([52|) 

In appendix-A, we work out the limiting behavior of Pn{x) as x ^ and x — > 1/A^. For general A^, one finds 



Piv(x) 



^r{N)r{Ny2) 

2JV-ir2(7v/2)r((Ar-l)/2) 
^jvA-^/2(l-A^x)(^'+^-4)/ 



-1/2 



as X 







as X 



1/N 



(55) 



Comparing this limiting behavior in the real case in Eq. ([55|l with that of the complex case in Eq. ([29|l one finds 
that while in the former Pn{x) diverges as a;""'^/^ as x — > 0, in the latter it approaches a constant. In the other limit 
X -^ l/N, both the densities approach zero as a power law (1 — Nx)" , but with different exponents ly = N'^ — 2 (for 
the complex case) and v — [N'^ + N — 4)/2 for the real case. 

100 n ' ^ ' ^ 1 



80 



complex case 
real case 



60 r> 



D. 




FIG. 1: The p.d.f Pn{x) of the minimum eigenvalue Amin vs. x for N — 4, for the complex and the real cases (Eqs. (|28p and 
(|52[) respecively) . In the complex case, the density approaches a constant as 2; — > 0, whereas for the real case, it diverges as 
x~^'^ as 2: ^ 0. 



Moments of Amin: One can use the explicit result for the p.d.f. Pn{x) of Amin in Eq. (|52p to calculate its fc-th 
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moment 

'■^/^ - --- .-- -- ■ ^- /iV + 2 iV-1 N'^ + N-2 l-Nx 



2 ' 2 ' 2 

r2 



dx 



/.i/iV / 

where we m.ade a change of variable y = — A^ + \/x in the first hne. We next use the following known integral [15| 

r x'<-'{x + z)-'^2F^{a,P,^,-x)dx^ ^^^^l\''.^^''f^'^~^^''K 
Jo r{a)r(a + fJ - ^ + a) 

in Eq. (|56p and also the value of A^r from Eq. (|53p to arrive at an explicit expression for the fc-th moment (valid for 
ah N), 

^^ ' r(iV/2)r(A: + iV2/2)r(fc+(7V + 3)/2)2^-i il^^'^/'^2' / ^^ 

One can verify that iJ-o{N) = 1, thus ensuring the correct normalization. For the average value of Amin we use k = 1 
and get 

Thus the expression for (Amin) for arbitrary N in the real case is considerably more complicated than its counterpart 
in Eq. pTjl for the complex case. One finds, from Eq. ([55)1 . that Hi{N) decreases with increasing N, e.g., /^i(l) = 1, 
^1(2) = (4— vr)/8, ^1(3) = (2— ■\/3)/9 etc. In appendix-B, we show that asymptotically for large N, ^i(iV) decays as 



^i{N)fii—; where c = 2 



Yerfc(l/V2) 



= 0.688641 • • • (60) 



where erfc(a;) = -?= /^ e "^ du is the complementary error function. The large N result in Eq. (j60p for the real 

case should be compared to that of the complex case where )J,i{N) — 1/N^. One sees that the average value of the 
minimum eigenvalue in the former case is less by a constant factor c = 0.688641 • • • compared to the later case. In 
Fig. 2, we plot both the exact formula for fii{N) in Eq. ((59)l and the asymptotic form in Eq. ((60|) against iV. 

V. CONCLUSION 

In this paper we have computed the exact probability distribution of the minimum eigenvalue Amin of an entangled 
random (both real and complex) pure state of a bipartite system composed of two subsystems whose respective 
Hilbert spaces have equal dimensions M = N. We have also computed exactly all the moments of Amin for all N. As 
a byproduct, we prove that (Amin) — l/N^ for all N for complex matrices, a result recently conjectured jI3|. The pdf 
of the minimum eigenvalue in the real case differs significantly from its complex counterpart. 

Apart from providing important informations on the nature of the entanglement of a random pure state as well as on 
the degree to which the dimension of the Hilbert space of a subsystem can be reduced, our result for the distribution 
of the minimum eigenvalue has some relevance in the general context of extreme value statistics. This subject has 
been around for a long time [23|, but has seen a recent resurgence due to its many applications in diverse areas 
such as engineering, economics and physical sciences [2g|. If the underlying random variables are independent and 
identically distributed then there are three possible limiting universal distributions for the extreme events, the Frechet, 
the Gumbel and the WeibuU distributions. However, much less is known when the underlying random variables are 
strongly correlated. In such cases, the limiting distribution (for large N) of the maximum is known exactly only in very 
few cases. For example, the limiting distribution of the largest eigenvalue oi a N x N Gaussian unitary random matrix 
(GUE) is given by the celebrated Tracy- Widom law [29|, which has found many recent applications [30] ■ Similarly, 
the Tracy- Widom law also describes the limiting distribution of the largest eigenvalue of Wishart matrices [3l|, Is^] , 
for random matrices with certain non-Gaussian entries [33| and the scaled height of a (1 -I- l)-dimensional growth 
models [30, [SJ] . The probabilities of large deviations of Amax, outside the regime of the validity of the Tracy- Widom 
law, have also been computed recently both for Gaussian [S^] and Wishart matrices [l^ [3l| . Other examples for 
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asymptotic c/N'' 




0.8 
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ln[N] 
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FIG. 2: A log-log plot of the exact formula of fii{N) in Eq. ((59} vs. A*' compared with the asymptotic formula /ii(A) ~ c/N^ 
vs. N with c = 0.688641 for the real case. 



vifhich the limiting distribution is known exactly include the maximum relative height of a class of one dimensional 
fluctuating interfaces in their steady states in a finite system [3a, [33] and 1//" noise signals i38]. In contrast, much 
less is known about the distribution of the extreme eigenvalues for finite N, a notable exception being the minimum 
eigenvalue for N x N Wishart matrices whose distribution was computed exactly by Edelman for all N 25]. In 
our present context, the eigenvalues of a random pure state are also strongly correlated due to the presence of the 
Vandermonde term Hixfc \^j ~ ^k\^ i'^ the jpdf ([T]). So our results provide another rare exactly solvable case for the 
distribution of the minimum of a set of N strongly correlated random variables, and this is not just for large N but 
for any finite N. 

Computing the distribution of Amin for unequal dimensions {M ^ N) of the Hilbert spaces of the subsystems 
remains a challenging open problem. 



APPENDIX A: LIMITING BEHAVIOR OF Pn{x) FOR THE REAL CASE 



In this appendix, we derive the behavior of P/v(a;) in Eq. (j55p . starting from the exact expression of PNi^) in 
Eqs. ((52)) and ([53)1 . The behavior near the upper limit x — > 1/A^ is simple to derive. Using 2-F'i(a,/3, 7, 0) = 1, one 
immediately finds from Eq. ([52ll that as a; — > 1/A^ 



Pn{x) ^ An N-'^'^ (1 - Arx)(^'+^-4)/2 



(Al) 



In contrast, deriving the behavior of Pn{x) as a; —+ is slightly more tricky. To derive this, we first use the following 
identity of the hypergeometric function [15] 



2Fi(a,/3,7,^) = (l-^)-^2i^i(/3,7- a, 7,^/(^-1)) 



to rewrite 



P^(^) = ^ (1 - iVx)(^'+^-4)/2 [1 - (iV - l)x]-(^-i)/2 2F1 ( 
VX \ 



N-lN^-4:N^+N-2 1- Nx 



2 ' 2 ' 
Now, in this form, it is easy to take the limit x -^ 0. One gets, as x — + 0, 



' 1- (iV-l)a; 



(A2) 



(A3) 



(A4) 
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Using further the foUowing identity [If 



2Fi(a,/3,7,l) 



r(7)r(7-a-/3) 



r(7-a)r(7-/3) 

and the expression for An in Eq. ([55)) we get, as a; ^ 

/^r(Ar)r(7VV2) 



Pw(a;) 



2W-ir2(Ar/2)r((iV-l)/2) 



'1/2 



(A5) 



(A6) 



APPENDIX B: ASYMPTOTIC BEHAVIOR OF iii{N) FOR LARGE iV FOR THE REAL CASE 

In this appendix we derive the asymptotic behavior for large N of /xi(iV) for the real case given in Eq. ([59)1 . We 
first use the following integral representation of the hypergeometric function [15J 



2Fi{a,P,j,z) 



1 



BiP,l-P) Jo 



\7-/3-l 



t^'-' {1 - ty-^'-' (1 - tzy^ dt, 



(Bl) 



where B{x,y) — T{x)T{ii)/T{x + y) is the standard Beta function. Using this representation, we can express fJ.i{N) 
in Eq. (|5^ as an integral 



"■("' = l^rii) r '"'('-""" I' + ("-W 



■dt. 



(B2) 



which is still exact for all N. Next we consider the integral above, rescale t — x/N and then take the large N limit 
as follows. 



t^^l - t)^/^ [1 + {N - l)t]'^ dt 



1 



m/2 



N 



xl/2 (1 - a;/iV)^/2 



N - 1 
l-\ X 

N 



dx 



oo ^l/2g-a;/2 



N^^^ Jo (1 + 2;)^ 



dx. 



(B3) 



Also, by Stirling's formula, r(A^) « \/2tt N^ ^/^ e ^ for large iV. Using these results in Eq. (jB2[) we get, to leading 
order for large N, 



MiW-]^ 



where the prefactor c is given by the expression 

4V2 r x^/^e^^^^ 



TT Jo (i + xy 



dx ^2 



1-Wyerfc(l/V2) 



= 0.688641 . . . 



(B4) 



(B5) 



where erfc(x) — -y= J e " dw is the complementary error function. 
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